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o 

fvj , Abstract. We prove Euler-Lagrange and natural boundary necessary opti- 

mality conditions for problems of the calculus of variations which are given by 
a composition of nabla integrals on an arbitrary time scale. As an applica- 
tion, we get optimality conditions for the product and the quotient of nabla 
variational functionals. 






1. Introduction 



u 
o 

^ ' A time scale T is any nonempty closed subset of K. The theory of dynamic 

+-J i Euler-Lagrange equations is a recent field attracting considerable attention — see 

^; [51[Sl[71[HJ[Il[THl[Tll[2Il[12]and references therein. For applications of the calculus 

of variations on time scales to economics see [TJ |3l [IT]- In particular, one obtains 

the well-known continuous [23 , discrete I15j, and quantum calculus of variations 

[5 by choosing T = M, T = Z, and q^" := {q''\k e No}, g > 1, respectively. 

^ ' This paper is dedicated to the study of general (non-classical) problems of cal- 

O^ . cuius of variations on an arbitrary time scale T. As a particular case, when T = M 

^— ' one gets the generalized calculus of variations |il2 with functionals of the form 

^. h( f f{t,x{t),x'{t))dt] 

f^ I where / has n components and H has n independent variables. Problems of calculus 

of variations of this form appear in practical applications of economics but cannot 
be solved using the classical theory (see [12] and the references therein). 
, , In the literature of the calculus of variations on time scales, the problems are 

rS I formulated and the results are proved in terms of the delta or the nabla calculus, 

j^ ■ Here we use a different approach. We make use of the duality technique recently 

introduced by Caputo |llj , obtaining the results for the nabla variational problems 
directly from the results on the delta calculus of variations. Such duality theory 
has shown recently to be very useful in control theory [24] . 

In contrast with |20| . we adopt here a backward perspective, which has proved 
useful, and sometimes more natural and preferable, with respect to applications 
in economics [TJ [2l |3]. The advantage of the backward approach here promoted 
becomes evident when one considers that the time scales analysis has important 
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implications for numerical analysts, who often prefer backward differences rather 
than forward differences to handle their computations. This is due to practical 
implementation reasons and better stability properties of implicit discretizations 

MM- 

The paper is organized as follows. In Section [2] some preliminaries on the recent 
duality theory on time scales are presented. Our results are then given in Section |3l 
we formulate the general (non-classical) problem of calculus of variations p.l|) on an 
arbitrary time scale; we obtain a general formula for the Euler-Lagrange equations 
and natural boundary conditions (Theorem 13.21) ; and interesting corollaries are 
presented for the product fCorollarv l3.4p and the quotient fCorollarv l3.6|) . Finally, 
in Section H] we illustrate the results of the paper with two examples. 

2. Preliminaries 

In this section we review some facts from |llj which we need for the proof of 
our main result (Theorem I3.2p . We begin by briefly recalling the basic definitions, 
notations, and facts concerning the delta and nabla differential calculus on time 
scales, which can be found in the monographs [Sl ITUlfTB] . 

A time scale T is an arbitrary nonempty closed subset of M. For each time scale 
T the following operators are used: the forward jump operator cr : T — > T, defined 
by a{t) := inf{s £ T : s > i} for t 7^ supT and cr(supT) ~ supT if supT < +00; 
the backward jump operator p : T — ;> T, defined by p{t) :— sup{s G T : s < i} 
for t ^ inf T and p(inf T) = inf T if inf T > —00; the forward graininess function 
/i : T ^ [0,00), defined by ii{t) :— a{t) — t; and the backward graininess function 
i^ : T — > [0,00), defined by ^{t) = i — p{t). A point i £ T is called right-dense, 
right- scattered, left-dense or left-scattered if a{t) — t, a{t) > t, p(t) — t, and 
p{t) < t, respectively. We say that t is isolated if p{t) < t < a{t), that t is dense 
if p{t) =t = cr{t). If supT is finite and left-scattered, we define T" := T \ {supT}, 
otherwise ¥":=¥. We say that a function / : T ^ R is delta differentiable at i G T** 
if there exists a number f^ (t) such that for all e > there is a neighborhood U of 
t such that \f{a{t)) - f{s) - f^{t){a{t) - s)\ < e\a{t) - s\ for aU s £ C/. We caU 
/^ (t) the delta derivative of / at t and / is said delta differentiable on T'' provided 
f^{t) exists for all t £ T". For f : T ^ X, where X is an arbitrary set, we define 
f^ := foa. A function / : T — > K is called rd-continuous if it is continuous at right- 
dense points in T and its left-sided limits exist (finite) at left-dense points in T. We 
denote the set of all rd-continuous functions by C^^ = Crd = Crd(T) = Crd(T;K). 
The set of functions / : T ^> E that are delta differentiable and whose delta 
derivative is rd-continuous is denoted by C^^ = C}^iT) = Cj^jj(T;R). 

In order to introduce the definition of nabla derivative, one defines a new set T^ 
which is derived from T as follows: if T has a right-scattered minimum m, then 
T„ = T\ {m\; otherwise, T^ = T. In order to simplify expressions, and similarly as 
done with composition with a, we define /''(i) := f{p{t)). We say that a function 
/ : T ^ M is nabla differentiable at t G T^ if there is a number /^ (i) such that for all 
£ > there exists a neighborhood U oit such that \fP{t)~f{s)-f^{t){p{t)-s)\ < 
e\p{t) — s\ for all s £ U. We call f^{t) the nabla derivative of / at t. Moreover, 
we say that / is nabla differentiable on T provided /^(i) exists for all t G T^. 
Let T be a time scale, / : T — >■ R. We say that function / is Id-continuous if 
it is continuous at left-dense points and its right-sided limits exist (finite) at all 
right-dense points. The set of all Id-continuous functions / : T — >• M is denoted by 
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Cjj — Cid ~ Cid (T) = Cid (T, M) , and the set of all nabla differentiable functions 
with Id-continuous derivative by C\^ — C\^ (T, M) . For the delta differential/integral 
calculus we refer the reader to ^ ; for the nabla differential/integral calculus we refer 
the reader to [13 Chap. 3]. 

Given a time scale T we define the dual time scale T* by T* := {s G R| — s G T}. 
If p and a denote the backward jump operator and the forward jump operator 
associated to T, then we denote by p and a the jump operators associated with T*. 
If V and /i denote respectively the backward graininess function and the forward 
graininess function associated to T, then we denote by v and p, the graininess 
functions associated to T*. 

Given a function / : T ^ M we define the dual function /* : T* — > M by 
/*(s) := /(— s) for all s e T*. For a given quintuple (T, p, cr, i/, /i) its dual will be 
(T*, /}, (7, j>, p) where p, ct, v^ p are given as follows: p(s) = —a{—s)^ a[s) = —p{—s), 
j>(s) — p*{s), and p(s) — h'*{s). 

Let / : T ^ K. The following holds: 

(i) If / is delta (resp. nabla) differentiable at to e T'' (resp. to G T^), then 
/* is nabla (resp. delta) differentiable at -to G (T*)^ (resp. -to G (T*)'') 
and /^(to) = -(/*)^(-to) (resp. f^ {to) = -(/*)^(-to), where A and V 
denote the delta and nabla derivatives for the time scale T; and A and V 
denote the delta and nabla derivatives for the time scale T* . 

(ii) / belongs to C^^ (resp. Cf^) if and only if /* belongs to C/^ (resp. C^^). 

Let a, 6 G T with a < b. We define the closed interval [a,b] in T by [a,b] :— {t G 
T : a < t < b}. Along this work we always assume that [a, b] denote an interval in 
a given time scale T. 

Let / : [a, 6] — >■ R be a rd-continuous (resp. Id-continuous). Then, 

/ /(t)At= / "r(s)Vs (resp. / f{t)Vt^ f "r(s)As) . 

Ja J-b \ Ja J-b J 

Lemma 2.1 ([H]). For a given Lagrangian L : [a,b]i^ x M. x M. ^ R the following 
identity holds: 



L{t,x''it),x^it))Vt = J \*(^s,{x*fis),{x*fis)) 



As 



for all functions X G Cl^{[a,b]), where the dual Lagrangian L* : [— 6, — a]"^ xRxR ^ 
M is defined by L*{s, x, v) = L{—s, x, ~v) for all (s, x, v) G [— fe, —a]^ x M x M. 

A note on the use of the term duality is in order. Duality theory is a standard 
concept in optimization theory that one learns in the earliest of courses beginning 
with its applications in linear programming and then in nonlinear programming. 
This theory has been extended to many settings, including optimal control theory 
and the calculus of variations. Typically, when one considers minimizing a function 
/(■) over a set C, the dual problem is one of maximizing a related function. In the 
present context, the word "duality" has a different meaning. The term duality does 
not refer here to the classical concept, well known to researchers in optimization, 
but to the one recently introduced in [11] (see also [24]). 
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3. Main Results 



Throughout we consider A,BeT with A < B. Now let [a, b] with a,b eT,b < B 
and a > A, be a subinterval of [A, S]. The general (non-classical) problem of the 
calculus of variations on time scales under our consideration consists of extremizing 
(i.e., minimizing or maximizing) 



b nb 



C[x]=HU .h{t,xP{t),x^{t))Vt,...,J Ut,xPit),x^{t))Vt\ 

{x{a) = Xa) (x(6) = Xb) 

over all x £ C/^. Using parentheses around the end-point conditions means that 
these conditions may or may not be present. We assume that: 

(i) the function H : R" — >■ K has continuous partial derivatives with respect to 
its arguments and we denote them by H^, i = 1, . . . ,n; 

(ii) functions (i, y, v) —>■ fi{t, y, v) from [a, 5] xR^ to M, z = 1, . . . , n, have partial 
continuous derivatives with respect to y, v for all t e [a, b] and we denote 
them by fiy, /„; 

(iii) fi, i = 1, . . . , n, and their partial derivatives are Id-continuous in t for all 

A function x G C^j^ is said to be an admissible function provided that it satisfies 
the end-points conditions (if any is given). The following norm in Cl^ is considered: 

||x||i= sup \x''{t)\+ sup \x^{t)\. 

tela.b] tela.b] 

Definition 3.1. An admissible function x is said to be a weak local minimizer 
(resp. weak local maximizer) for (|3.ip if there exists S > such that C[x] < C[x] 
(resp. C[x] > C[x]) for all admissible x with ||a; — x\\i < 6. 

The next theorem gives necessary optimality conditions for problem p.ip . 

Theorem 3.2. If x is a weak local solution of problem p.ip . then the Euler- 
Lagrange equation 

n 

(3.2) ^i/a-^i[5;],...,-^n[i])(/.^(t,5''W,i^W)-/.«(t,i''W,i^(t))) = o 

holds for all t e [a, 5]^, where Ti[x] — J fi{t, x''{t),x'^ {t))Vt, i ~ 1, . . . ,n. More- 
over, if x{a) is not specified, then 

n / a(a) 

(3.3) ^ij;(J-i[i],...,J-„[5;]) / f,y{t,5:P{t),i^(t))Vt 

-/™(a(a),x''(a(a)),i^(a(a)))j =0; 
if x{b) is not specified, then 

n 

J2H^{^l[i],■■■,^n[i])f^v{b,i''{b),i^{b)) = 0. 
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Proof. Since i is a weak local extremizer for p.l|) . it follows by Lemma uA\ that x* 
is a weak local extremizer for the dual problem 



-a 

„*\A/ 



C*[x*] =Hl^j ^ /r(s, (x*)-(s), (a;*)^(s))As, 

{x*{-a) = Xa) {x*{-b)^Xb) 

over all x* € C'rdd"^' ~^])- Applying [201 Theorem 3.2] we conclude that x* 
satisfies the following conditions: 

n 

j2 Hii^iii*],- ■ ■,Km) {{f:.f{^, {i*ns), {i*f{s)) 

i=l 

for all s e [-6, -a]«, where J'*[x*] = /J^" f*{s, {S:*f{s), {i*)^{s))As, i = l,...,n. 
Moreover, if x*{—b) is not specified, then 

n 

(3.4) Y. H[{:Fl[i% . . .,K[i*])IU-b, {i*f{-b), {i*)^{-h)) = 0; 
1=1 

if x*{—a) is not specified, then 

(3.5) J2 H[{:FI[x% . . .,K[i*]) ( /»(/5(~a), {iT{p(-a)), {i*f{p{-a))) 



By Lemma I^TTl 



Jp{-a) J 



•^rri- / /*(s,(F)'^(s),(r)-(s))As 

-b 
l-b 



i = 1, . . . , n, and t G [a, &]«. From the duality of the Euler-Lagrange equations |lll 
Theorem 6.10] it follows that 

ifL)^ (^, (s^Ti^), (i*)^(^)) - fty (^, (s^Tis), (^ )^(^)) 

= fZ {t,i''{t),i^{t)) ~ Uy {t,if'{t),x^{t)) , 

i = 1, . . . ,n. This establishes relation p.2p . Now assume that a;(6) is not specified. 
By duality, condition (13. 4p holds. Since 

/*,(-&, (x*)n-6),(r)'^(-&)) = -/„(6,i^(6),i^ (6)), 
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i = 1, . . . ,n, we have 

n 

^ HiiJ^.ii], . . . , J-„[S])/„(&, iP{b),i''{b)) = 0. 

If x{a) is not specified, then by duality condition p.Sp holds. Observe that 

/*,(p(-a),(rf(p(-a)),(r)^(p(-a)) = -/„(a(a),i''(a(a)),i^(a(a)), 
and 

J p{ — a) Ja 

From the above it follows p.3p . D 

Choosing T = K in Theorem 13.21 we immediately obtain the following result: 
Corollary 3.3. (Th. 3.1 and Eq. (4-1) in [ll]^ If x is a solution of the problem 

C[x]^h( fi{t,x{t),x'{t))dt,..., fnit,x{t),x'{t))dt] — > extr 

{x{a) = Xa) (x(6) = Xb) 
then the Euler- Lagrange equation 

5^F,'(.Fi[i],...,J-„[5]) /,,(i,x(i),i'(t))-^/™(t,£(i),i'(i))) -0 

holds for all t G [a, b], where J-i[x] — J fi{t, x{t), x'(t))dt, i — 1, . . . , n. Moreover, 
if x{a) is not specified, then 

n 

^H[{Ti[x\,...,J^n[x])ftv{a,x{a),x'{a)) = 0; 

i=l 

if x{b) is not specified, then 

n 

Y.H'AJ'iii], ■ ■ ■,Tn[i])fUb,m,i'{b)) = 0. 
Corollary 3.4. If x is a solution of the problem 

-CN = ( / fiit,xP{t),x^it))yt] I f f2it,xP{t),x^{t))yt] -^ extr 



{x{a) = Xa) ix{b) = Xb) 
then the Euler- Lagrange equation 

T2[i] {fl{t,iP{t),i''{t))-f,y{t,iP{t),S:''{t))) 



+ Fi[x] {fZit,xP{t),i^it))~f2yit,xPit),S:^{t)))^0 
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holds for all t G [a, 6]„. Moreover, if x{a) is not specified, then 

/ Ma) \ 

J'2[i] U fiy{t,i''it),x^{t))\7t~fUaia).S:'{aia)),x^iaia)))\ 

( /-^(a) \ 

+ -^i[£]iy /2,(i,x''(t),£^(t))Vi-/2.(a(a),f''(a(a)),i^(a(a)))j = 0; 

if x(b) is not specified, then 

J'2[x]fUb,xP{b),x^ (b)) + Ti[x]f2v{b,i''{b),S:^ (b)) = 0. 



Remark 3.5. In the particular case T = M, Corollary 13.41 gives a result of [T2]: the 
Euler-Lagrange equation associated with the product functional 



"1^' ^ (/ 



flit,xit),x'it))dt] I f f2it,xit),x'it))dt 



IS 



J^2[x] ifly{t,x{t),x'{t)) - —fi^{t,x{t),x'{t)) 

+ ^lN (f2y{t,x{t),x'{t))^j^f2.{t,x{t),x'{t))\ =0 

and the natural condition aX t = a, when x{a) is free, becomes 

^2[x]fiv{a,x{a),x'{a)) + Fi[x]f2v{a,x{a),x' {a)) = 0. 

Corollary 3.6. If x is a solution of the problem 

J^f,it,xPit),x^it))Vt 
jl\x\ = -^ > extr 

!lf2{t,xp{t),x^{t))m 

{x{a) = Xa) {x{b) = Xb) 
then the Euler-Lagrange equation 

fl{t,iP{t),i''{t))--f,y{t,iP{t),S:^{t)) 

- Q {fl{t,iP{t),i^{t)) - f2y{t,iP{t),i^{t))) = 

holds for all t £ [a,b]^, where Q = j} ^ . Moreover, if x{a) is not specified, then 
fiy{t, x''{t),i^{t))Vt - A, {a{a), iP{a{a)), i^ {<j{a))) 

r(a) \ 

f2y{t, iP{t), i^{t))Vt - /2, {a{a),iP{a{a)), i^ {a{a))) J = 0; 

if x{b) is not specified, then 

/i,(5, S:P{b),x^{b)) - Qf2y{b, xP{b),x^{b)) = 0. 
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Remark 3.7. In the particular situation T = M, Corollary 13.61 gives the following 
result of [12j : the Euler-Lagrange equation associated with the quotient functional 

^^^j ^ Slh{t.x{t),x'{t))dt 
jlh{t,x{t),x'{t))dt 
is 

hy{t,x{t),x'{t))-Qf2y{t,x{t),x'{t)) 

d 



^^ ,{h,{t,x{t),x'{t)) - Qf2y{t,x{t),x'm = 

and the natural condition sd, t ~ a, when x{a) is free, becomes 
/ii,(a,x(a),x'(a)) - Qf2via,x{a),x' (a)) = 0. 

4. Examples 
Example 4.1. Consider the problem 

a;(-l) = 1, 2;(0)=0. 
If i is a local minimum of ()4.ip . then the Euler-Lagrange equation must hold, i.e., 

(4.2) 25;^^(i)Q2 + Qi = 0, 
where 

/O pO 

{x^{t)fVt, Q2^F2[i\-- tx^{t)Vt. 

If Q2 = 0, then also Qi = 0. This contradicts the fact that on any time scale a 
global minimizer for the problem 

Fi[x\ = / {x^{t)fVt — > min 

a;(-l) = 1, a;(0)=0 

is x{t) — ~t and !Fi[x\ = 1. Hence, (52 7^ and equation (|4.2p implies that the 
extremals for problem (|4.ip are those satisfying the delta differential equation 

(4.3) iVV(,)^_|i 

subject to boundary conditions x{—l) = 1 and x(0) = 0. A solution of ()4.3p 
depends on the time scale. Let us solve, for example, this equation on T = R and 
onT= {-l,-i,0}. On T = R we obtain 

Substituting (|4.4I) into functionals J-"! and 7^2 gives 

( 48Ql+Ql _ ^ 
(4.5) <^ 48Q^ '^l 

L 24Q2 ^2- 
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Solving the system of equations (|4.5p we obtain 





Therefore, 

x{t) ^ -f - 2t 

is an extremal for problem (|4.1I) on T = M. 
The solution of g^ on T = {-1, -5,0} is 

if t = -1 
(4.6) ^(t) = <;i+^ ift=-i 

if t = 0. 




Constants Qi and Q2 are determined by substituting (j4.6p into functionals J-i and 
J-2 . The resulting system of equations is 

(4.7) 

Since system of equations (|4.7|) has no real solutions, we conclude that there exists 
no extremizer for problem (|4.ip on T = { — 1,-^,0} among the set of functions 
that we consider to be admissible. 

Example 4.2. Consider now the problem 

/°,(xV(i))2Vt 



^N = TIT 



(4.8) /_2(a;^W + (a^^W)')Vt 

x(-2)=4, x(0)=0. 

If a; is a local minimizer for (|4.8|) . then the Euler-Lagrange equation must hold, i.e., 

= [2x^(t) - Q(l + 2i^(t))]^, t e [-2, 0]„, 



where 



Therefore, 






As x{-2) = 4 and a;(0) = 0, we have Q ^ 1. Thus i'^^ (t) = 0. The solution of the 
delta differential equation x'^^ (t) — 0, x{—2) = 4, a;(0) = 0, does not depend on 
the time scale: x{t) ~ —2t is an extremal for problem (j4.8p . 
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